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Abstract 

The problem of describing optical turbulent structures, arising in resonant media with 
high Fresnel numbers, is reviewed. The consideration is based on the probabilistic ap- 
proach to pattern selection, ascribing a probability distribution of patterns to systems 
with multiple possible structures. The most probable structure corresponds to the mini- 
mal expansion exponent, or to the minimal expansion rate, characterizing the phase-space 
expansion of a dynamical system. Turbulent photon filamentation is studied. The most 
probable filament radius and the number of filaments are found, being in good agreement 
with experiment. 
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1. Experiments on Optical Structures 

Spatio-temporal optical structures appear when electromagnetic waves propagate through 
a medium possessing some sort of nonlinearity. One can distinguish three types of such media, 
depending on the physical origin of nonlinearity. 

(i) Kerr medium. This is a passive medium with a nonlinear relation between polarization 
and electric field, when susceptibility depends on the field intensity. Many liquid crystals 
pertain to this kind of matter. Different patterns arising in the Kerr medium are reviewed in 
Refs. [1-3]. 

(ii) Photorefractive medium. This pertains to a medium with a strong photorefractive 
effect. The photorefractive effect is the change in refractive index of a medium resulting from 
the optically induced redistribution of electrons and holes. The modulation of the medium 
refractive index is usually realized as a four-wave mixing process. Two input optical fields 
produce intensity grating, shifting electrons from donor atoms, which induces a space charge 
wave, equivalent to the formation of a ripple field. The nonlinear interaction of the input fields 
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with the medium leads to the generation of two output waves. An example of a photorefractive 
medium is the bismuth selicon oxide crystal Bii2Si02o- Various patterns are reviewed in Ref. 
[2,3]. 

(iii) Resonant medium. This, evidently, is an active medium typical of laser systems. Nonlin- 
earity arises because of an effective interaction of resonant atoms through the common radiation 
field. Different spatio-temporal structures were observed in a Na2 laser, CO2 laser, dye laser 
and in many semiconductor and vapour lasers, as is reviewed in Refs. [3-6]. 

Pattern formation in nonlinear optics, which is also called optical morphogenesis, has many 
similarities for different media. This especially concerns the case when nonlinear optical samples 
have the same geometry. The most often used is the cylindrical geometry typical of lasers. 
Characteristic structures of cyhndrical samples are filaments, aligned along the cyhndrical axis, 
where the radiation intensity is much higher than in the space between these filaments. The 
appearance of such structures is easily noticeable in the transverse cross-section, where the 
high-intensity filaments are exhibited as bright spots. For brevity, we call these filamentary 
structures photon filaments [5,6]. 

In the case of cylindrical samples, the peculiar properties of photon filaments essentially 
depend on the Presnel number 



in which R is the cylinder radius, L is the characteristic length, and A is the radiation wave- 
length. The Presnel number is a sort of an aspect ratio comparing the transverse and longitu- 
dinal characteristics. 

The Presnel number in optics plays the same role as the Reynolds number in hydrodynamics. 
When increasing the Reynolds number, the fluid passes from laminar motion to turbulent one. 
A similar transition happens in optics when increasing the Presnel number, and occurs around 
F ^ 10. 

At small Presnel numbers F < 5, there can exist several transverse modes with a regular 
behaviour. These transverse structures are regular in space, forming ordered geometric arrays 
seen as polygons in the transverse cross-section, with the number of bright spots proportional 
to F^. The structures are also regular in time being either stationary or periodically oscillating. 
The type of such regular structures is prescribed by the sample geometry and corresponds to 
the empty-cavity Gauss- Laguerre modes imposed by the cylindrical geometry. This regular 
behaviour in optics is analogous to the laminar motion in hydrodynamics, and it is well un- 
derstood theoretically both for nonlinear media, such as a Kerr medium and photorefractive 
crystals [1-3], and for laser media [4]. 

Por large Presnel numbers F > 15, the appearing spatio-temporal patterns are principally 
different from the empty-cavity modes. The modal expansion is no longer valid and the bound- 
ary conditions have no importance. The medium looks as a bunch of independently flashing 
fllaments, whose number is proportional to F. The fllaments are chaotically distributed in the 
transverse cross-section, are not correlated with each other, and are flashing apcriodically in 
time. This spatio-temporal optical chaos is similar to hydrodynamic turbulence, because of 
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which it is called optical turbulence. Since this phenomenon is characterized by the formation 
of bright filaments with a high photon density, it has been named [7] turbulent photon filamen- 
tation. This kind of optical turbulence was observed in nonlinear media, such as photorefractive 
crystals [2,3] as well as in different lasers, such as CO2, dye, and vapour lasers (see discussion 
in Refs. [5-7]). 

The theory of turbulent photon filamentation was advanced in Ref . [7] . This theory is based 
on the probabilistic approach to pattern selection [6] . 

2. Probabilistic Pattern Selection 

The problem of pattern selection arises when evolution equations possess several solutions 
corresponding to different spatio-temporal structures. This problem is rather common for 
many nonlinear systems, optical media being only a particular case [8]. Therefore, let us, first, 
formulate a general approach to pattern selection, relevant for any dynamical system. 

Let us consider a set of functions, in general complex, yk{x,t) enumerated by the index 
k = 1,2, . . . and depending on a collection of variables x G D C M'^ and on time t G M. The set 
of = {yk{x,t)} is the dynamical state. For the compactness of notation, it is convenient to 
denote by one index i the pair k and x, writing y{t) = {yi{t)}, with yi{t) = yk{x,t). Then the 
dynamical state y{t) can be treated as a column with respect to i. In this notation, the system 
of evolution equations can be presented as 

j^y{t) = viy,t), (1) 

with a velocity filed being also a column v{y,t) = {vi{y,t)}. Evolution equations containing 
only the first derivatives in time are called to have the normal form. Any system of equations 
containing higher time derivatives can always be reduced to the normal form [9]. 

The dynamical state y{t) C pertains to a phase space T — ®iTi, being a tensor product 
of the spaces Ti D yiit). An elementary phase volume of the phase space T can be written [6,7] 
as \bV{t)\, with 

bV{t)^Y^bym. (2) 

i 

The temporal behaviour of the phase volume element (2) is characterized by the expansion 
exponent 

5V{t) 



a{t) = In 



(3) 



which shows whether the phase volume expands with time or contracts according to the law 

\6T{t)\ = |5r(0)|e"W . 

The expansion exponent (3) can be presented in several equivalent forms [6,7]. Introducing 
the multiplier matrix M{t) = [Mij{t)] with the elements 

M«(t) . MO , (4) 

6yj{0) 
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we can transform Eq. (3) to 

a{t) ^]n\det M{t)\ . (5) 

Another form of the expansion exponent (3) can be obtained by means of the expansion 
matrix X{t) — [Xij{t)] with the elements 

^'^^^^ = -J^ ■ 

The variation of the evolution equation (1) yields 

I M{t) = X{t)M{t) . (7) 

^Prom Eqs. (6) and (7), we find 

ait) = Re /* TxXit') dt' . (8) 

^0 

If the multiplier matrix (4) possesses eigenvalues Hn{t) enumerated by a multi- index n, then 
Eq. (5) gives 

a(t) = $:in|/i„(i)| . (9) 

n 

Let us introduce the local expansion rate 

A{t) ^ i a{t) . (10) 

The latter, depending on the form of the expansion exponent, can be written in the following 
three ways: 

A{t) = - In detM(t) , (11) 
i 

A{t) = -Re rTrX(t') dt' , (12) 



t Jo 

m = T.7 i^Kwi ■ (13) 

^Prom Eq. (13), we have 

limA(t) = ^A„, (14) 

n 

where the right-hand side is the sum of the Lyapunov exponents 

An = lim J In |Af„(i)| , (15) 
provided the corresponding limit exists. 



It is easy to check that 

lim|.W = X:A„, (16) 

which immediately follows from Eqs. (9) and (15). Since from Eq. (8) one has 

a{t) = Re TrX{t) , 

then, according to Eq. (16), one gets 

hm Re TrX(t) = V . 

t— >oo 

n 

These limiting properties remind us those of the Gibbs entropy 

Sait) ^- j p{y, t) lnp{y, t) dy , 

in which p{y, t) is a probability density. For the Gibbs entropy in steady states, one has [10-14] 
the limit ^ 

The quantity dSo/dt is called the entropy production rate or, simply, the entropy production. 

Therefore, the expansion exponent (7{t) is a kind of dynamic entropy variation, and the expan- 
sion rate A(t) resembles the entropy production rate. More details on the entropy production 
can be found in Refs. [11-16]. 

Notice that normally an invariant probability measure of a steady state is singular [13] and 
therefore the Gibbs entropy So{t) is Sait) — — oo, which is physically unacceptable. 

It is worth mentioning that in dynamical theory the quantity Re TrX(t) is sometimes termed 
the contraction rate. According to Eq. (12), the expansion rate (10) could also be named the 
contraction rate. The usage of the words "expansion" or "contraction" is rather philological. 
But, probably, it is more logical to call A{t) the expansion rate, since A{t) > means really 
the expansion of the phase volume |5r(i)|, while A{t) < imphes the contraction of the latter. 

Now let us turn to the problem of pattern selection, when the evolution equation (1) pos- 
sesses a set of solutions corresponding to different spatio-temporal structures. If we label these 
different patterns by an index a, then for each pattern we have a solution y{a,t) and, respec- 
tively, an expansion exponent a{a,t) and an expansion rate A{a,t). The problem is how to 
classify these different solutions and the related patterns. 

The probabilistic approach to pattern selection [6,7] is based on the idea that each pat- 
tern can be characterized by its probabilistic weight. That is, there should exist a pattern 
distribution p{a,t). The latter can be derived by minimizing the pattern information 

Ip{t) = J p{a,t)lap{a,t) da + J p{a,t)a{a,t)da , (17) 
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consisting of the sum of the Shanon information and lost information, under the normahzation 
condition 

j p{a, t) da = 1 . 

This imphes the minimization of the conditional pattern information 



ip{t) = Ip{t) + l{t) j p{a, t) da - 



in which l{t) = ha.Z{t) — 1 is a Lagrange multiplier. Minimizing Ip{t) results in the pattern 
distribution ^ 

P{a,t) = ex.p{-a(a,t)} , (18) 

with the normalization factor 

Z{t) — J exp{—a{a,t)} da . 
Taking into account the relation (10), we obtain 

p(«, ^) = ^ {-^(«' ^) ^} ■ (19) 

In this way, the most probable pattern corresponds to the minimal expansion exponent a{a,t) 
or, equivalently, to the minimal expansion rate A{a,t), that is, one has the relation 

minA(Q!,t) <^=^ msixp{a,t) . (20) 

Thus, we come to the following conclusion. 

Principle of Pattern Selection. The most probable pattern for a dynamical system at 
a given time is defined by the minimal expansion rate. 

Let us emphasize that the formulated here approach of Probabilistic Pattern Selection is a 
very different from the Prigogine [17-21] principle of minimal entropy production in station- 
ary states. The latter principle deals with the entropy production S = dS/dt defined as time 
derivative of thermodynamical entropy. It says that the entropy production S{a) in a station- 
ary state is minimal with respect to internal thermodynamic variables a. The basic points, 
distinguishing the Probabilistic Pattern Selection from the Prigogine principle, are as follows: 

(i) This approach is developed for arbitrary dynamical systems. In particular, it may con- 
cern the evolution equations for thermodynamic variables or for statistical averages of some 
operators. 

(ii) The theory is valid for any dynamical states, not necessarily stationary ones. 

(iii) Not solely the most probable pattern is defined, but the probability distribution for all 
feasible patterns is suggested. 
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3. Turbulent Photon Filament at ion 

Now we shall apply the theory of Section 2 to the problem of turbulent photon filamentation, 
discussed in Section 1. Consider a system of N resonant atoms interacting with electromagnetic 
field as is described by the Hamiltonian 

H^Ha + Hf + Haf (21) 

consisting of the terms: The Hamiltonian of two-level resonant atoms 



Ha = E^o[^ + Stj , (22) 

where loq is a transition frequency and is a spin-1/2 operator. The radiation-field Hamilto- 
nian 

with electric field E and magnetic field H = V x A. The vector potential is assumed to satisfy 
the Coulomb gauge cahbration V • A = 0. The atom-field interaction Hamiltonian 

^ /I \ 
^a/ = -E -J*-Ai + D,-EoJ , (24) 



i=l 



in which Aj = A(r, t), and the transition-current and transition-dipole operators are 

Ji = lujo {dS+ - d*Sr) , Bi = dS+ + d*Sr , 

where d is the transition dipole, d = doe^; Sf are the rising or lowering operators, respectively; 
and Eoi is a seed field. 

We shall consider the evolution equations for the functions 

u{r, t) = 2< S-{r, t) > , s(r, t) = 2 < S'{r, t) > , (25) 

being the statistical averages of the corresponding operators. Introduce the notation 

f{r,t)^fo{r,t) + fra4M (26) 

for an effective field acting on an atom, with the cavity seed field 

/o(r,i) = -2id-Eo(r,i) (27) 

and the radiation field 



3 p 

frad{r, t) = - - i-fp J 



jfco|r-r'| -ifeo|r-r'| 

n - 4 ^- u*{r',t) 

/co r — r Ko r — r 



dr' , (28) 
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where ko = coo/c, p is the density of atoms, and 7 = AkQdl/3 is the natural width. The seed 
field 

Eo = - El e*('=^-'^*) + - El e-*^'^^-'^*) (29) 
selects the longitudinal mode with the frequency cu — kc satisfying the quasiresonance condition 

|A| 

— <C 1 , A = u} — u}o. 

The evolution equations for the atomic variables (25) can be derived [5] by eliminating the 
field variables and invoking the semiclassical approximation. This yields 

— = -{iuo + -f2)u + fs , = -272|m|2 + {u*f + f*u)s , 

(9s 1 

Q-^=-^{u*f + ru)-^,{s-C), (30) 

where 71 is the level width, 72 is the line width, and > is a pumping parameter. 

Without the loss of generality, we may look for the solutions to Eqs. (30) in the form of a 
bunch of Nf filaments, 

Nf Nf 

u{r, t) = J2 Mr±, t)e''' , s{r, t) = J2 ^n(r±, t) , (31) 

n=l n=l 

with r± = y/x'^ + y'^. In the case of very small Presnel numbers, F <C 1, there is the sole 
filament, so that Nf — 1 and it„ and Sn are uniform in space. For low Presnel numbers F < 10, 
when the regime of laminar filamentation occurs, an expansion over Gauss-Laguerre modes is 
more appropriate [4]. And for large Fresnel number F > 10, when the regime of turbulent 
filamentation develops, there appear many uncorrelated filaments. Then the functions Un and 
Sn are essentially nonzero only around the axis of an n-th filament, but fastly decrease outside 
the filament, so that UmUn ~ Smn, SmSn ~ Smn, and UmSn ~ Smn- Here we are interested in the 
latter case of turbulent photon filamentation. 

To describe the turbulent regime, let us introduce the averaged functions 

Unit) = ^ j Un{r±,t) dr , s„(t) = / Sn{r±,t) dr , (32) 

with the averaging over a cylinder enveloping the n-th filament. The volume of the enveloping 
cylinder is \4 = vra^L, with being the cylinder radius and L, the sample length. The radius 
of the enveloping cylinder On is related to the filament radius r„ by the conservation-energy 
equation 

j |«n(r±,t)|' dr = Vn\un{rn,t)\^ . (33) 



8 



If the filament profile is well approximated by the normal law exp(—rj_/2r^), the relation (33) 
gives 

an = 1.82 rn • (34) 

The sample is assumed to have the cylindrical shape of radius R and length L. For these 
values and for the radiation wavelength A, the inequalities 

typical of lasers, are valid. For the functions (32), characterizing the turbulent regime, we 
obtain the equations 

ds 1 

= -gnl2\Un\^ - 2 «/l + fl'^n) " ll{s„ - C) , (35) 

where 

Tn = 72(1 - Qns) , Q„ = a;o + 72^^, , /i = -id ■ Eic"''"* , 
and the efective coupling parameters are 

37P f sm{kor-kz) , f cos{kor - kz) 

9n = -;— r dr , g^= — dr . 36) 

472 Jvn kor 472 Jvn kor 

The equations for different filaments are naturally decoupled since the filaments are not corre- 
lated with each other. 

The nonlinear system of equations (35) can be solved by means of the scale separation 
approach [5,22]. The details of this solution can be found in Refs. [6,7]. Calculating the 
expansion rate (12) and minimizing it with respect to the filament radius, we obtain the most 
probable radius 

r/ = 0.3v^ , (37) 

which is in good agreement with experiment, as is discussed in refs. [5-7]. The formula (37) 
has been specially checked in the regime of turbulent photon filamentation of CO2 and dye 
lasers and found to well describe the experimental observations [23-27]. 

The most probable number of filaments can be estimated from the normalization integral 



^Jsdr^C, (38) 



where the integration is over the whole volume of the sample V = ttR^L. Considering the 
developed turbulent regime, when the population difference inside each filament of radius Vf 
has reached the value of the pumping parameter and the population difference outside the 
filaments is Sgut, we have from the integral (38) the equation 

NfVfC + {V- NfVf)sout = CV . 
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^Prom here, we find 



With the filament radius (37), this gives 




(39) 



Nf^AF , ^ ~ , (40) 

which also is in good agreement with experiment. 

Filamentation in laser media is a rather general process and should arise for different sam- 
ples, provided the Presnel number is sufficiently high. This should concern as well the so-called 
random lasers [28]. 

The theory presented above has to do with the well developed optical tm4:)ulcncc corre- 
sponding to large Fresnel numbers. An interesting question is how to describe the intermediate 
region, where 5 < F < 15, and the arising turbulence is yet weak, being intermittent with the 
remnants of the regular filaments, characterized by the Gauss-Laguerre modes. Then between 
the laminar and turbulent filaments there could appear a kind of Josephson interference, as it 
happens for electromagnetic vortices [29,30]. In the frame of the probabilistic pattern selection, 
described in Sec. 2, we could treat laminar and turbulent filaments as independent, calculating 
their probabilities. In the intermediate region of Presnel numbers, where these two types of 
filaments coexist, their probabilities should be comparable. At low Fresnel numbers, the prob- 
ability of the laminar filament must prevail, while at high F, the probability of the turbulent 
filaments is to become much larger. In this way, the whole picture for arbitrary F, varying in 
the interval [0, oo\ can be done in the frame of the probabilistic pattern selection. 
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